Chaotic dot-superconductor analog of the Hanbury Brown Twiss effect 
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As an electrical analog of the optical Hanbury Brown Twiss effect, we study current cross- 
correlations in a chaotic quantum dot-superconductor junction. One superconducting and two 
normal reservoirs are connected via point contacts to a chaotic quantum dot. For a wide range 
of contact widths and transparencies, we find large positive current correlations. The positive cor- 
relations are generally enhanced by normal backscattering in the contacts. Moreover, for normal 
backscattering in the contacts, the positive correlations survive when suppressing the proximity ef- 
fect in the dot with a weak magnetic field. 



In quantum theory, identical particles are indistin- 
guishable. Under exchange of any pair of particles, the 
many-body wavefunction remains invariant up to a sign, 
positive for bosons and negative for fermions. This ex- 
change symmetry was used in the pioneering interferom- 
eter experiment with photons, by Hanbury Brown and 
Twiss . Several theoretical works have suggested 
different analogs of this experiment with electrons in 
mesoscopic multiterminal conductors. It has been shown 
that the fermionic statistics of electrons leads to neg- 
ative correlations between currents flowing in different 
terminals. Such negative correlations were also recently 
observed experimentally |^]. 

When normal conductors are connected to a super- 
conductor, correlations are introduced between electrons 
and holes due to Andreev reflections at the normal- 
superconductor interface, a phenomenon known as the 
proximity effect. The influence of the proximity effect 
on the current auto-correlations, i.e. the shotnoise, in a 
two-terminal diffusive normal-superconductor junctions 
was recently studied HJ^. 

In multiterminal conductors, Andreev reflection can 
lead to positive cross correlations between currents flow- 
ing in the contacts to the normal reservoirs |^,^. So far, 
positive correlations have been predicted only for single 
mode junctions ^-p^. Moreover, in multiterminal dif- 
fusive junctions it was found that cross correlations are 
negative in the absence of the proximity effect [ pl| . 

This raises two important questions: i) are the positive 
correlations in normal-superconducting junctions a large 
effect, of the order of the number of modes in multimode 
junctions, and if this is the case ii) is the proximity effect 
necessary to obtain these positive correlations? In this 
paper we give an answer to these two questions. The posi- 
tive correlations are large, and surprisingly, get enhanced 
by normal backscattering at the normal-superconducting 
interface. Moreover, positive correlations can exist even 
in the absence of the proximity effect, if the normal- 
superconductor interface is nonideal. 

We study the current correlations in a system consist- 
ing of a chaotic quantum dot connected via point con- 
tacts to one superconducting and two normal reservoirs. 



Systems consisting of chaotic dots coupled to supercon- 
ductors have recently attracted a lot of interest. The 
generic properties of the model makes our result qualita- 
tively relevant for multiterminal normal-superconducting 
structures with random scattering. 
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FIG. 1. A chaotic quantum dot (grey shaded), acting as a 
beam-splitter, is connected to two normal reservoirs (A^'i and 
A'^2) and one superconducting reservoir {S) via quantum point 
contacts. 

A schematic picture of the system is shown in Fig. |^. 
A quantum dot is connected to two normal reservoirs 
(iVi and N2) and one superconducting reservoir {S) via 
quantum point contacts. The contacts to the normal 
and superconducting reservoirs have mode independent 
transparency Tn and Ts respectively and support TV and 
M transverse modes. The conductances of the point 
contacts are much larger than the conductance quanta 
2e^//i, i.e. iVFjv, MFs > 1, so Coulomb blockade effects 
in the dot can be neglected. The two normal reservoirs 
are held at the same potential V and the potential of the 
superconducting reservoir is zero. 

We consider the case where the classical motion in 
the dot is chaotic on time scales longer than the er- 
godic time Terg- The quasiparticle dwell time in the 
dot, %l Etvi, is assumed to be much larger than Tg^g, but 
much smaller than the inelastic scattering time. Here 
Eth = (SiVFjv + MTs)S/-K, where S is the mean level 
spacing in the dot. Under these conditions random ma- 
trix theory can be used to describe the statistical 
properties of the scattering matrix Sd of the dot. The 
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scattering matrix S of the combined system, dot and su- 
perconductor, can be expressed [Q in terms of the scat- 
tering matrix Sd of the dot and the Andreev reflection 
amphtude at the contact-superconductor interface. 

Due to the random scattering in the dot, the cur- 
rent Ii(t) in contact i fluctuates around its quantum 
statistical average 7^. We study the zero- frequency 
spectral density of the current cross-correlations P12 — 
2jdt A/i(i)A/2(0) where AIj(t) = U{t)-h. The corre- 
lation can be expressed in terms of the scattering matrix 
S. We consider the limit of zero temperature and a po- 
tential eV much lower than Exh and A, where the energy 
dependence of the scattering matrix S can be neglected. 
Moreover, at \E\ < A no quasiparticles can escape into 
the superconductor and the scattering matrix S describes 
only scattering between the normal reservoirs. 
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where S°j are matrix amplitudes {N x N) for injected 
quasiparticles (e or h) of type (3 in lead j to be back 
reflected as quasiparticles of type a in lead i. 

Noting that the current fluctuation is just the sum of 
the fluctuations of electron and hole currents, the noise 
power can be conveniently written S , 
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where the noise power for correlation between quasi- 
particle currents, is given by 



Pri=aV- 
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[S^tyS^j^iS^^^yS^t] , (3) 



with a = +(— ) for a ^ (3 (a ^ f3). The expression 
for the correlations in Eqs. is an extension of 

the result |^ for a purely normal conductor, taking into 
account both electron and hole quasiparticles. However, 
unlike the normal cross correlations, which are manifestly 
negative, the cross-correlation P12, can be positive, be- 
cause the correlations between different types of quasi- 
particles, Pf2 + Pyz^ ^'^^ positive. Note however that 
P12 + P12 - {P12 + P12), the cross-correlator between 
probability currents, is manifestly negative a conse- 
quence of the ferminonic statistics of the quasiparticles. 

The ensemble averaged correlations (P12) are calcu- 
lated using the statistical properties of the scattering ma- 
trix Sd of the dot. We first consider the case with ideal 
contacts rjv,r5 = 1 and no magnetic field in the dot, 
i.e. conserved time reversal symmetry. In this case it is 
useful to decompose Sd as (see e.g. Ref. Q) 
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where U{U') is a unitary matrix of dimension 2N x 
2N {AI X M), uniformly distributed in the unitary 



group. The diagonal matrices r2N,2N and Tm.m have 
min(2iV, M) elements y/1 — r„ and the rest unity {t2N,M 
and tM,2N follow from the unitarity of Sd)- Here T„ are 
the transmission eigenvalues, which have a density p^ ] 
p[T) = {2N + M)/(2^)VP - r™„/(T\/r^), where 
T,mn = [{2N - M)/{2N + M)f is the smallest possible 
eigenvalue. 

Inserting the decomposition into the expression for 
the total scattering matrix S we find Q the quasi- 
particle scattering amplitudes S"^ in Eq. (|^) as e.g 
S"=e = Uy/T^/{2 - T)U* (T = 1 - r2^,2Ar)- Summing 
over injection contacts in the expression for the individual 
quasiparticle current correlators in Eq. (|^) and inserting 
the scattering amplitudes we get 

P?l ^ v^ixWil - 7^)^7*ClC/^7^C/^C2] (5) 
h 

and similar for the other terms P12 ■ Here TZ — T^/ (2 — 
T)^ are eigenvalues of the matrix product {S'^'^Y S'^^ and 
the matrix Ci is diagonal with elements {Ci)n = 1 for 
n < N and otherwise. The matrix C2 = 1 — Ci. The 
ensemble average of Eq. (||) is carried out in two steps. 
First, P"^ is averaged over the unitary matrix U, using 
the diagrammatic technique in Ref. ||l5|. This gives to 
leading order in N, M (i.e neglecting weak localization 
corrections) for the total current corellations in Eq. (||), 

(P„),,.Ki£!(t,|K(i-R,i-^22|^).(e) 

We then perform the average over transmission eigenval- 
ues by integrating {Pi2)u weighted by the density p{T). 
Using that to leading order in N we have (tr(7^)tr(l — 
7^))T = (tr(7^))T(tr(l-7^))T, we get (< .. >=< .. >c/,t) 
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where /(7) = 1 + 67 + 7^ and 7 = 2N/M, the ratio be- 
tween the conductances of the point contacts connected 
to the normal and the superconducting reservoirs. The 
correlation is normalized with Pq = AVNe^/h, i.e. it is 
large, of order N. The correlation (P12), plotted in Fig. 
H, is positive for a dominating coupling to the supercon- 
ductor, but crosses over (at 7 « 0.5) to negative values 
upon increasing the coupling to the normal reservoirs. 

In the general case, it is difficult to provide a simple 
explanation for the sign and magnitude of the correla- 
tions (P12), since they result from a competition between 
{Pi2) + {Pi2') and {Pi2) + {Pii), which have opposite sign 
and generally are of the same magnitude. It is however 
possible to get a qualitative picture in certain limits by 
studying the different contributing scattering processes. 
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FIG. 2. The ensemble averaged current cross-correlation 
(P12) with (solid) and without (dashed) proximity effect in 
the dot, as a function of 7 = 2N/M. The correlations in the 
presence of the proximity effect are positive for a dominating 
coupling to the superconductor (7 < 0.5), but are negative 
for all 7 when suppressing the proximity effect with a weak 
magnetic field in the dot. 

In the limit 7 — > 0, the couphng to the normal reser- 
voirs is negligible and a gap in the spectrum opens up 
around Fermi energy in the dot As a consequence, 
quasiparticles injected from one normal reservoir are An- 
dreev reflected, efi^ectively direct at the contact-dot inter- 
face with unity probability back to the same reser- 
voir. There is thus no partition of incoming quasiparticles 
and hence no noise. Pi 2 = 0. 

Increasing the coupling to the normal reservoir, the 
probability of normal reflection(~ 7) as well as cross An- 
dreev reflection {'-^ 7^), from one reservoir to the other, 
becomes finite. Since normal reflection is the dominant 
process, we can neglect the terms in P12 in Eq. (||) con- 
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taining the cross Andreev reflection amplitude (e.g. S 
This gives (P12) = W{e^ /h){tY{Sll{Sll)^ Sfl{Sl'l)^)) = 
2Po7, positive since only terms in (Pf^) -I- (Pi''2'^) con- 
tribute. This expression shows that the correlations 
can be explained as partition noise of injected electrons, 
which have probability ^ 1 to be Andreev reflected (effec- 
tively at the contact-dot interface) and probability ~ 7 
to be normally reflected. 

In the opposite limit, 7 3> 1, the coupling to the super- 
conductor is weak and the Andreev reflection probability 
is small, of order I/7. The correlation can be written as 
(P12) = ^lW{e^/h){iT{Sll{Sll)^Sl^{St,n = -P0/7, 
which is negative since we find that in this limit only 
contribute. This shows that the 



terms in (P 
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correlation can be explained as partition noise of elec- 
trons which have probability ^ 1/2 to be normally re- 
flected (without reaching the dot-superconductor con- 
tact) and probability ~ I/7 to be Andreev reflected. 

When the proximity effect is suppressed by a weak 
magnetic field in the dot, it is no longer possible 
to express the correlations directly in terms of the trans- 



mission eigenvalues T, as in Eq. We instead use 

the fact that the scattering matrix Sd of the dot itself is 
uniformly distributed in the unitary group, without con- 
straints on symmetry of Sd- Since the scattering matrix 
amphtudes S^'P in Eq. (|) can be expressed in terms 
of Sd, the correlations P^2 can be directly averaged over 
Sd using the diagrammatic technique in Ref. This 
gives for the total correlations 
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plotted in Fig. |. The gap in the spectrum is suppressed 
and the scattering processes responsible for the positive 
correlations in the presence of the proximity effect are 
strongly modified, most notably in the limit 7 <IC 1. As a 
consequence, the correlations are manifestly negative for 
all 7, i.e. there are no positive correlations in the absence 
of the proximity effect, similar to what was found for a 
metallic diffusive system 11 1. 

Until now we have only considered ideal point contacts, 
with FatjE^ = 1. In an experimental situation, it is 
often difficult to obtain an ideal contact between the dot 
and the superconductor. It is therefore of interest to 
study the situation with a nonideal dot-superconductor 
contact. We consider first the case with proximity effect 
in the dot. To calculate the current correlations in this 
case we note that a nonideal interface changes the 
density of transmission eigenvalues, p(P), but not the 
distributions of the unitary matrix U in Eq. (^). The 
transmission eigenvalue density is calculated numerically 
]l6| , p^ for different contact transparencies Fg < 1 and 
the integrals in Eq. (^) are then performed. 

The resulting correlations are plotted in Fig. ||. Sur- 
prisingly, the main effect of normal backscattering at the 
dot-superconductor contact is to cause a crossover from 
negative to positive correlation for a strong coupling to 
the normal reservoirs. In this limit, 7^1, injected 
quasiparticles undergo at the most one scattering event at 
the dot-superconductor contact before leaving the junc- 
tion. The expression for the scattering matrices S'^J^ in 
Eq. (^) simplify considerably and we can derive an ana- 
lytical expression for the correlations, giving |17 
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where Reh = F|/(2 — Fs)^ is the Andreev reflec- 
tion probability of quasiparticles incident in the dot- 
superconductor contact. There is a crossover from neg- 
ative to positive correlations already for Pe/t = 1/2, i.e 
F5 — 2(V2 — 1) w 0.83, in agreement with the full nu- 
merics in Fig. ^. 

Since F5 < 1, now R^h is smaller than one. As a conse- 
quence, there is additional partition due to the possibility 
of normal reflection at the dot-superconductor contact. 
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The partition noise discussed above, from electrons be- 
ing either Andreev reflected or normally reflected (with- 
out reaching the dot-superconductor contact) , is thus re- 
duced to —PaRehll- However, the additional normal 
reflection at the dot-superconductor contact give rise to 
a noise term 2Po^e/i(l — Reh)/l, with opposite sign [to- 
gether they give Eq. (||)]. For Rf-h < 1/2, the second 
term is dominating, causing the crossover to positive cor- 
relations. 
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FIG. 3. Ensemble averaged current cross correlations 
Pi2 in the presence of normal backscattering at the nor- 
mal-superconducting interface (see inset). The contact trans- 
parencies Fs are 1 (solid), 0.8 (dash-dotted) 0.6 (dashed) and 
0.4 (dotted). The correlations are plotted as a function of 
2N/{Mrs)- Note that the correlations for a dominating cou- 
pling to the normal reservoirs crosses over from negative to 
positive on increasing the normal backscattering. 

Interestingly, in the absence of a proximity effect in 
the dot, we find in the same way that the correlation 
for 7 1 is also given by Eq. (||), the argument being 
the same as in the presence of the proximity effect. The 
proximity effect thus plays no role in this limit, where the 
quasiparticles undergo at most one Andreev reflection. 
This shows that the proximity effect is not a necessary 
condition for positive correlations. 

Finally, we note that the effect of normal backscatter- 
ing in the contacts between the dot and the normal reser- 
voirs is to enhance the positive correlations for a domi- 
nating coupling to the superconducting reservoir. In the 
limit with tunnel barriers in all contacts, the transmis- 
sion eigenvalue density is known analytically and the 
correlation (P12) follows from Eq. (p|) p8t| 
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where 7 — 2A^FAr/(MFs). The correlations are thus pos- 
itive for 7 < (\/5— 1)/2 and 7 > {\/5 + l)/2, and negative 
for intermediate values. 

In conclusion, we have studied the current cross- 
correlations in a three terminal superconducting-chaotic 



dot analog of the Hanbury Brown Twiss interferometer. 
We find that the correlations are positive for a wide range 
of junction parameters, and survive even in the absence 
of a proximity effect in the dot. The magnitude of the 
positive correlations is large, proportional to the num- 
ber of transport modes in the contacts to the dot, which 
should simplify an experimental observation. 

We acknowledge discussions with C.W.J. Beenakker, 
W. Belzig, E. Sukhorukov, and G. Johansson. The work 
was supported by the Swiss National Science founda- 
tion and the program for Materials with Novel Electronic 
Properties. In the final stages of this work we became 
aware of Boerlin et al. , who investigate a three ter- 
minal junction in the timnel limit. In this limit, our 
results coincide. 
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